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Abstract. In Higgs-otic inflation a complex neutral scalar combination of the hP and 
MSSM Higgs fields plays the role of inflaton in a chaotic fashion. The potential is protected 
from large trans-Planckian corrections at large inflaton if the system is embedded in string 
theory so that the Higgs fields parametrize a D-brane position. The inflaton potential is then 
given by a DBI-I-CS D-brane action yielding an approximate linear behaviour at large field. 
The inflaton scalar potential is a 2-field model with specihe non-canonical kinetic terms. 
Previous computations of the cosmological parameters (i.e. scalar and tensor perturbations) 
did not take into account the full 2-field character of the model, ignoring in particnlar the 
presence of isocurvature perturbations and their coupling to the adiabatic modes. It is well 
known that for generic 2-field potentials such effects may significantly alter the observational 
signatures of a given model. We perform a full analysis of adiabatic and isocurvature per¬ 
turbations in the Higgs-otic 2-field model. We show that the predictivity of the model is 
increased compared to the adiabatic approximation. Isocurvature perturbations moderately 
feed back into adiabatic fluctuations. However, the isocurvature component is exponentially 
damped by the end of inflation. The tensor to scalar ratio varies in a region r = 0.08 — 0.12, 
consistent with combined Planck/BICEP results. 
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1 Introduction 

It is very hard to protect fundamental scalars in field theory from large quantum corrections. 
It is well known that these quantum corrections easily drive the scalar masses to the size of 
the largest ultraviolet cut-off in the theory. Even so, light scalars seem to play an important 
role in our understanding of the observed properties of particle physics and cosmology. On 
the one hand the recent observation of the Higgs particle at the LHC [1] is consistent with a 
fundamental SU{2) x U{1) complex doublet remaining light well below any ultraviolet cut¬ 
off. On the cosmology side, precision data increasingly favours the existence of an inflaton 
scalar, with mass well below the Planck scale, leading to a period of fast expansion of the 
universe during which quantum perturbations generate CMB anisotropies and plant the seeds 
for galaxy formation. 

Low-energy supersymmetry provides a rationale for the stability of the Higgs mass 
against quantum fluctuations. However, the fact that no trace of SUSY has been found yet 
at the LHC and the measured value of the Higgs mass rrih ~ 125 GeV cast some doubts 
on the idea of low-energy SUSY. Indeed, such a Higgs mass is quite high compared to the 
expectations from the MSSM and can only be reached for very heavy SUSY spectra. One 
can think of disposing altogether with the idea of SUSY and admit that the Higgs field is 
fine-tuned to be light, perhaps on the basis of anthropic arguments. Still this is probably not 
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enough, since if we insist on the validity of a non-SUSY SM structure at higher energies, the 
renormalisation group evolution of the Higgs self-coupling is such that the potential becomes 
unstable (or metastable) for scales larger than ~ 10^^ GeV, possibly indicating the presence 
of new physics at these scales, well below the Planck scale (see e.g.[2]). 

A most elegant way to stabilise the Higgs potential for scales above ~ 10^® GeV is again 
SUSY. Indeed, the scalar potential of e.g. the MSSM is positive definite for scales above 
SUSY breaking and hence the potential is automatically stable. It was shown in [3, 4] (see 
also [5-7]) that under very general assumptions, if SUSY is broken at scales GeV, 

the resulting Higgs mass at low-energies is consistent with its measured mass ~ 125 GeV. 
Of course, with such a high SUSY breaking scale the Higgs mass is again unprotected and 
one should admit fine-tuning as the ultimate reason for its lightness. This fine-tuning could 
come from an underlying multiverse within string theory, which would also be at the root of 
the understanding of the smallness of the cosmological constant. 

On the cosmology side, a fundamental scalar degree of freedom, the inflaton, seems to 
provide the simplest explanation for a variety of cosmological observations. That scalar again 
will have to be protected from acquiring a large mass and interactions that spoil slow-roll. 
On the other hand, as we already remarked, light scalars usually do not stay light in a field 
theory. Intuitively, one may expect having two fundamental scalars light to be even more 
unlikely than a single (Higgs) one. Thus if the Higgs field itself could act as the inflaton the 
combination Higgs/inflaton would look more likely within a landscape of theories. Of course, 
independently from any landscape argument, it is a natural question whether the Higgs and 
inflaton scalar could be one and the same. This has been considered often in the past, see 
e.g. [8-15], see [16] for a review and references. 

In ref. [17] it was proposed that the neutral Higgs system of the MSSM with SUSY 
broken at a large scale ~ 10^^ GeV could be in charge of cosmic inflation. This goes under 
the name of Higgs-otic inflation [18], since it is in some sense a MSSM Higgs version of 
Linde’s chaotic inflation [19]. This proposal is quite economical since it addresses several 
issues simultaneously. It provides stability for the Higgs scalar potential at the right scale, 
is consistent with the observed value of the Higgs mass and a neutral Higgs component 
acts as a complex inflaton. The inflaton has a trans-Planckian field range and leads to 
a certain variety of (non-canonical) 2-field chaotic like inflation. In order for this to be 
sensible we need a theory of quantum gravity, which is naturally identified with string theory. 
The Higgs/inflaton vev corresponds to the position of a D-brane moving on periodic cycles 
in a version of monodromy inflation [20-25], [26],[27-37], [38-40]. The potential for the 
Higgs/inflaton is dictated by the Dirac-Born-Infeld -|- Chern-Simon action of the moving D- 
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brane. Fluxes induce a scalar quadratic potential for the Higgs/inflaton system, with masses 
of order the SUSY breaking scale Mss- However, the kinetic term of the scalars is not 
canonical but rather is determined by the scalar potential itself. This leads to a flattening 
of the scalar potential at large fields [41], with a leading linear behaviour at large inflaton. 
In general, a substantial contribution of tensor perturbations is generated. 

The Higgs-otic model is a 2-field inflaton system. In [18] a study of the cosmological 
observables was made for the model focusing only on the adiabatic perturbations and ignoring 
possible 2-field specific effects like the generation of isocurvature perturbations. This was 
an important pending issue since 2-field effects can, in principle, substantially modify the 
cosmological observables and furthermore the Planck satellite has provided stringent bounds 
on isocurvature perturbations. In the present paper we perform a systematic analysis of the 
observables in the Higgs-otic 2-field inflation system. We find that, as expected, adiabatic 
and isocurvature perturbations form a coupled system and there is super-horizon evolution 
of the curvature perturbations. This leads in general to a relative increase of adiabatic 
perturbations and consequently to a reduction of the tensor to scalar ratio r compared to the 
computation in ref. [18]. The range of variation of Ug gets smaller and is centered around 
the region alowed by Planck data with a tensor to scalar ratio in a range r = 0.08 — 0.12. 
Moreover, the isocurvature component is always very suppressed at the end of inflation, 
consistent with upper Planck bounds. 

The structure of this paper is as follows. We review the main points of Higgs-otic 
inflation in the next section, in which the relevant definitions and the inflaton potential 
are described. In section 3 we review the main issues of 2-field inflation as applied to the 
case of Higgs-otic inflation. Section 4 presents the results for Higgs-otic inflation for three 
representative points in the parameter space of the induced soft-terms for the Higgs-inflaton 
system. The latter are determined by a real positive parameter 0 < ^ < 1, with A = 
{mjj — -|- m|), h being the neutral Higgs scalars driving inflation [18]. The first 

case {A = 0.83) correspond to the canonical Higgs-otic model in which the lightest scalar 
field at the minimum of the potential (at scale Mss) can be identified with the SM Higgs 
field. A second case {A = 0.7) analyses how those results are changed if there is some effect 
(like modified RG running) slightly modifying the Higgs-otic setting. For completeness we 
finally present a third case with A = 0.2 in which the inflaton cannot be identified with the 
MSSM Higgs fields but could be relevant in some extensions of the MSSM. At the end of this 
section we show the expectations for the r — Ug plot in the Higgs-otic model for a variety of 
initial conditions and mass scales. Finally, section 5 is left for the conclusions. 
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2 Higgs-otic Inflation 


In this section we briefly review and extend the results presented in [17, 18] on Higgs-otic 
inflation. We thus motivate the form of the non-minimal 2-fleld inflaton potential to be 
analysed in the following sections. Readers interested only in the inflationary analysis may 
safely jump to the next section. 

Higgs-otic inflation refers to theories in which the inflaton is a complex scalar giving 
rise to gauge symmetry breaking, while attaining large field inflation. The most obvious and 
natural candidate for that is the SM Higgs field itself, as described in [18] . Nevertheless the 
same idea may be applied to other BSM fields introduced for other purposes, as we briefly 
discuss below. The essential ingredient is the identification of a complex inflaton with the 
position moduli of some Dp-brane system in string compactifications. ^ The motion of the 
brane corresponds to the gauge symmetry breaking through a scalar vev. For concreteness 
we will consider the position moduli of a DT-brane describing cycles in a torus inside 
a Type HB orientifold compactiflcation, although the setting may be easily generalized to 
other string configurations As we emphasised above, the reason to go to a string setting 
in order to implement this idea is twofold. First, we will be interested in producing a 
large held model in which inflation implies trans-Planckian excursions of the inflaton. In 
order to do so in a sensible manner we need a theory of quantum gravity. Indeed within 
string theory the vevs of scalar fields may be trans-Planckian and one can still maintain an 
effective potential Lagrangian which makes sense. Second, string theory has the required 
properties/symmetries in order to keep under control Planck-supressed corrections which 
become obviously important at trans-Planckian held values. There are modular and shift 
symmetries which force these potential corrections to be subleading. 

So let us consider, for definiteness, a Type HB orientifold yielding a theory with N = 1 
supersymmetry in = 4. Gauge interactions appear from stacks of Z17-branes wrapping 
4-cycles in the CY compact manifold. Those will include the SM gauge interactions as well 
as the SM fields. They come from zero modes of = 8 complex scalars and A upon 
reduction to D = 4. In particular, the position of the Z17-branes is parametrized by the vevs 
of zero modes. On a stack of N D7 the gauge group is U{N) and the scalars transform as 
adjoint chiral multiplets. The gauge symmetry may be reduced to that of the SM (or some 
extension) by e.g. an orbifold twist in the compactiflcation. Generically only some fields of 
the adjoint scalars <I> survive, like those parametrising the motion of those D7-branes which 

^The same structure appears in terms of continuous complex Wilson lines in models with extra dimensions. 
From the string theory point of view they correspond to equivalent T-dual compactifications. 

^For an introduction to orientifold constructions see [42] and references therein. 
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can move along flat directions. This is the case of L)7-branes associated to SM Higgs vevs in 
ref. [18]. 

The generic presence of closed string fluxes in Type IIB compactifications will give rise to 
a non-trivial potential for these fields, see ref. [43-50]. We thus assume, as is customary, that 
there are imaginary self-dual (ISD) fluxes G 3 acting as a background. Such backgrounds are 
known to be solutions of Type IIB D = 10 equations of motion in warped CY backgrounds. In 
these backgrounds there are two classes of ISD fluxes, with tensor structure G'{o, 3 ) and <-^( 2 , 1 )) 
respectively. The first correspond locally to components Gjgg and the second to Gjgs- From 
now on we will denote G = G 123 and S = esjkG^j^. The first class of flux, G, breaks SUSY 
and gives rise to SUSY-breaking soft terms, scalar and gaugino masses. The second class, 
S, preserves SUSY and gives rise to supersymmetric F-terms, i.e., //-terms. These induced 
masses will eventually be identified with the mass scale of the inflaton. The cosmological 
bounds on the density scalar perturbations applied to our model fix this scale to be of order 
10^^ — 10^^ GeV [51]. In order to identify the inflaton with an MSSM Higgs boson we will see 
that the simultaneous presence of both classes of fluxes, G and S, is required. This implies 
that the supersymmetry breaking scale Mgs (ie. the scale of the soft terms) in this large field 
inflation setting will also be of order 10^^ — 10^^ GeV. This is consistent with the scenario 
of Intermediate SUSY breaking in [3, 5-7] which was shown to be consistent with a Higgs 
mass of 126 GeV [4]. The structure of mass scales in the Higgs-otic setting is summarised in 

fig-(l)- 
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Figure 1. Energy scales in the Higgs-otic Inflation scenario [18]. Below 10^^ GeV the light degrees 
of freedom in the Higgs sector are given by the SU{2) doublet Hl- Above this scale SU{2) is broken 
and they lie within the neutral components of h and H. 

The dynamics of the fields living on the Z17-branes is described by the Dirac-Born- 
Infeld (DBI) plus the Ghern-Simons (GS) actions [42]. Inserting the above two classes of flux 
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backgrounds G, S in the DBI+CS actions one obtains an effective D = A bosonic action [18] 

Cm = STr - VF(ci>) - , (2-1) 

V ^9ym ^ / 

where there is a non-canonical kinetic factor /($) which depends on the scalar (F-term) 

potential Vp as 

/(‘h) = 1 + , (2.2) 

with 

Z~‘^o 

Vp($) = Tr |G*^> - . (2.3) 

and Str denotes the symmetrised trace over gauge indices. Here gs is the string coupling, 
V 4 the volume wrapped by the F7-brane and Z~‘^ = 1 in the absence of warping, which 
we will assume in what follows. In addition there is the D7 tension gj = M®/(27r)’^, with 
Ms = the string scale. To obtain Eq. (2.1) we have kept all the terms involving 

higher order powers of the position modulus as required in a consistent large field inflation 
model. This is doable because the DBI action contains all perturbative corrections in a', 
whose effect is encoded in the non-canonical kinetic factor /(<k). For small field / ~ 1 and 
we recover the result prescribed by the 4d supergravity effective action [18]. 

The above action concerns the U{N) adjoint in the world volume of F7-branes. How¬ 
ever, it still applies after we make an orbifold projection that converts the adjoint into a set 
of bifundamentals, some of which are identified with the Higgs field. For example [17, 18] one 
can start from a set of 6 F7-branes yielding an initial U{6) gauge and project out yielding 
a U{3) X U{2) X F(l) gauge group. One linear combination of the 17(1)’s can be identified 
with hypercharge, whereas the other two are anomalous and become massive in the usual 
way. The adjoint contains doublets surviving projection, i.e. 


4 > = 


O3 

O2 Hu 

y Hd 0 J 


(2.4) 


where Hu and Hd can be identified with the usual MSSM Higgs fields. Plugging this decom¬ 
position in Eq. (2.3) and taking the trace one obtains the scalar potential ^ 

Vf = [([G] - \S\f\h\^ + ([G] + \S\)^\H\^] , (2.5) 

^The same scalar potential may be obtained in an = 1 supergravity formulation not resting on the 
DBI-bCS action. However, such a SUGRA treatment will miss the important effect of the non-canonical 
kinetic term, which gives rise to the flattening of the potential, see [18]. 
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where one defines 


h = 


72 


H = 


Vi 


( 2 . 6 ) 


with 7 = TT — Arg(G5') given by the relative phase of both fluxes. 

The model is essentially double chaotic inflation with the key difference of non-canonical 
kinetic terms. To see different limits of this expression it is useful to define the real variable 
A which controls the relative size of SUSY (S) versus non-SUSY (G) fluxes, i.e. 


2|5G| 

|G|2 + |5|2 • 


(2.7) 


Note that 0 < A < 1. This parameter may also written in terms of the masses of the above 
defined neutral scalars, 


mu 

mh 


1 + A 
1 - A 


( 2 . 8 ) 


There are two interesting limits in which this 2-field scalar potential becomes effectively a 
single field potential: 


• A = 1. In this case |G| = IS"! and the held h becomes massless. If we want to eventually 
hne-tune a Higgs doublet to remain at low-energies as a SM Higgs, much below the 
inhaton mass scale, we would need to be close to that situation. The Higgs held H will 
be the one producing inhation in this limit. 


• A = 0. This happens, in particular, in the SUSY case in which G = 0. In this case the 
inhaton mass would be supersymmetric and a = -|- \H\‘^ acts as an effective 

single inhaton. In this case one can hardly identify any of these scalars with (neutral) 
SM Higgs helds, since they will have large masses, of order the inhaton mass ~ 10^^ 
GeV. On the other hand such helds could be identihed with other BSM helds which 
were not required to survive at low energies and could act as inhatons. Examples of 
these could be SU{5) scalar triplets or SU{2)ji doublets of left-right symmetric models. 

The most economic option is of course identifying the helds with the MSSM Higgs 
doublets. In such a case we would need one full SM Higgs doublet to remain massless 
below the SUSY breaking scale ~ 10^^ GeV. This requires hne-tuning of the huxes G, S. 
However, one has to take into account the running of masses in between the UV scale (the 
compactihcation/string scale Ms) and the SUSY breaking scale Mgs — 10^^ GeV. Taking 
into account this effect it was found in ref(...) that one needs A ~ 0.83. Despite this fact, 
in our inhationary analysis we will present results for a range of possible A values. It should 
be kept in mind that only results in which A is not far from a central value A ~ 0.83 are 
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consistent with the identification of the infiaton with MSSM Higgs fields. Other values of A 
could correspond to other BSM alternatives, as discussed above. 

It is useful to write down the scalar potential in Eq. (2.5) in angular variables. Defining 
the D-term fiat direction as 


a = \Hu\ = \Hd\ , Hu = , (2.9) 

one rewrites the potential as 

H((t,0 ) = M|<j(^l-Hcos0)fj2 , (2.10) 

where we define 

M|5 = , |Gp = Z-2(I4^7)-i(|G|2 + |,S|2) , (2.11) 

and 0 = 9 — Arg{GS). Note that Mss measures the size of the soft terms and hence the 
infiaton mass. One can estimate the size of these parameters from purely stringy arguments 
by assuming an isotropic compactification, Z ~ ~ 1 and standard GUT values for the 

gauge couplings [18]. One finds 

|G| = [Z-V4-'/V7^ , (2.12) 

Mis = V4fi79s\Gl ~ 0.05Mf\Gl , (2.13) 

where Mg is the string scale, assumed to be of order 10^® — 10^^ GeV. In this way one has 
an infiaton mass of order Mss — 10^^ GeV, consistent with the cosmological bounds. 

In angular variables the relevant piece of the action may be written as 

CM = f{a,e) (2{d^af + - Mlsil - Acose)a^ . (2.14) 

with f{(J,9) = 1 + ^{V 4 g, 7 gs)~^V{a, 9). Eq.(2.10) will be our infiaton potential, although 
one should not forget that the kinetic terms are not canonical and this fact should be taken 
into account. This will lead to a general flattening effect, as we discuss below. Note that 
the shape of the scalar potential (before flattening) depends essentially on the value of the 
H-parameter defined above. In fig. 2 we show the structure of the scalar potential for three 
characteristic values A = 0.1,0.5,0.95. 

In the one-field limits, H = 0,1, one can always get canonical kinetic terms by making 
the following field redefinition 





(2.15) 
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Figure 2. Scalar potential for three different values of A, A = 0.1 (upper left), A = 0.5 (upper right) 
and A = 0.9 (below). 


In those limits one has a potential 

Vf = , (2.16) 


with (/) = (T for ^ = 0 (G = 0 or S' = 0) or (/> = \H\ for ^ = 1 (|G| 
is then 



From Eq. (2.15) one finds the canonically normalised scalar ip 


|S|). The kinetic factor 
(2.17) 


If = 


^|<^|^2 + | G | 2 | 0|2 + 


^|G|-^sinh-M|G||0|/^/2] . 


(2.18) 


There is not a simple analytic expression for the inverse function (p{ip), which is what one 
would want in order to get the final form of the canonically normalised potential. However, 
for large helds one can check that the second piece in Eq. (2.18) is small and one can invert 
analytically ip{4>) and obtain the scalar potential for the new variable 


Vf = {gsV^P7)\GfH^ ~ (<7sV4ft7) (^1 + 8|G|V - l) • (2-19) 

Note that, for large field, the potential becomes linear in the inflaton (p. This behaviour 
for large field is also present in the results of the complete 2-field analysis below. However, 
finding such an analytic expression is not always possible. In general, for multiple fields there 
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does not exist a transformation that makes the metric flat on the moduli space parametrized 
by (/>“. Therefore, to make a complete analysis we will have to use the general Lagrangian of 
multiple helds 

Ud = - V{ct>) , (2.20) 

where in our case the metric will be given by 

4/(a,0) 0 \ 

( 2 . 21 ) 

0 j 

and use the generalized expressions of the cosmological observables for non-canonical kinetic 
terms and multiple fields introduced in the next chapter. 

3 Basics of 2-field infiation 

We are interested in studying cosmological perturbations of a system of two scalar fields 
coupled to Einstein gravity. The action for such a system, if one allows for a curved field 
space, takes the form (SvrG = Mp^ = 1): 

s = J V^d^x - VicP)^ , ( 3 . 1 ) 

where is the spacetime metric with determinant g and Gab the held-space metric. The 
scalar fields (j)^ and ([P' span a 2-dimensional manifold with all relevant information contained 
in Gab- The Christoffel symbols, Riemann tensor and Riemann scalar of the scalar manifold 
are defined in the usual way. Assuming that the scalar fields depend only on (cosmic) time, 
())“ = </>o(t) and that the spacetime metric is the flat FRW metric: 

ds^ = —dP -|- a^{t)dx^dxi , (3.2) 

the equations of motion for the scalar fields are given by: 

'ci>l + = 0 . (3.3) 

This can expression can be written more economically using the covariant derivative defined 
as: = d(^Q The function a is the scale factor and the Hubble parameter, H, 

is defined as: H = c?dn(a). The equation of motion for the scale factor is: 

^' = ^^<^0 + ^) , (3.4) 

where = Gab<Po^o- Solutions to Eqs. (3.3) and (3.4) determine the background evolution 
of the system. 
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It is useful to define the slow-roll parameters: 

H 


rj = 


He 


(3.5) 

(3.6) 


Physically, e measures the change of the Hubble sphere of the local universe. Hence, 0 < e < 1 
corresponds to a shrinking Hubble sphere, which is a good definition of an inflationairy period. 
In order for inflation to last a sufficient amount of time the derivative of e has to be small 
as well, so r/ <C 1. It is possible to use the equations of motion (3.3) and (3.4) to relate e to 
(pQ = Gab^Pofpo, which implies that rj is related to the tangential acceleration of the background 
trajectory. 

In order to make connection with observations it is necessary to consider quantum per¬ 
turbations around classical background solutions. This involves expanding the scalar and 
metric degrees of freedom in terms of the background quantities {4 >q and g))^) and perturba¬ 
tions, finding the equations of motion for the gauge invariant perturbations and solving them 
[52]. The spectral index and tensor-to-scalar ratio are defined in terms of the power spectra 
of the quantum fluctuations of scalar and tensor modes as we will shortly review. 

It is convenient to define a local frame on the trajectory in field space given by [53-55] : 


= <^o) > 

90 

N = —,— • {~G22<Po ~ Gi24>0, Gh^q Gi2^q) ■ 


(3.7) 

(3.8) 


y/G<po 

These vectors form an orthonormal basis of the tangent space of the scalar manifold. Hence, 
they can be used to decompose the physically relevant quantities along the normal and 
tangential directions with respect to the background trajectory. In particular, the derivatives 
of the scalar potential can be written in this basis as 1^ = T^daV and Vjv = N^daV (o labels 
both components of the basis vectors as well as helds). The total acceleration of the fields 
can be derived by taking the covariant derivative of r“, using the equations of motion and 
projecting with A^“. One obtains: 


(po = -2>H(po -14 


(3.9) 


The parameter rj needs to be generalised to capture the full dynamics in field space. Recalling 
that in single field settings rj measures the acceleration of the scalar field, one finds that its 
convenient multifield generalisation is 

1 


V = 


HcPo 


Dcp- 


(3.10) 
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which measures the acceleration for the field 4>°‘. 

In the local (T, N) basis, 77 “ is projected on the following parallel and perpendicular 
components: 


V\\ 

v± 


(t>o 

Vn 


such that: 


(3.11) 

(3.12) 


= ^T^ + v±N^. 


(3.13) 


r/y measures the tangential acceleration, responsible for the variation of the modulus of the 
background trajectory velocity whereas rj± measures the normal acceleration causing the 
background trajectory to curve. 

The equations of motion for the perturbations are given in terms of the gauge-invariant 
Mukhanov-Sasaki variables [56, 57]: 

+ (3.14) 

which can also be decomposed using the basis of Eq. (3.8) as : 

= aTaQ^ , (3.15) 

= aNaQ^ . (3.16) 


Quantisation of these variables is done by the standard procedure, defining the canonical 
momenta and demanding the correct commutation relations. The variables are given in 
Fourier space by: 

= J + (3.17) 


where Oq, and ah are the usual creation and annihilation operators. We have switched to 
conformal time defined by: dr = l/a{t)dt. The Greek index a labels the quantum modes of 
the perturbations, consistency requires that a = 1,2 in the 2-field case. 

The equations of motion for the Mukhanov-Sasaki variables in this basis take the form 

[55]: 


'2„,T 


d'^v. 


dr'^ 

d^v^ 

dr"^ 


+ 2aHrij_ 


dv, 


N 


z tt‘2 2 T I 

-a H + 

2 tt2 2 N 

-a H - 


— 2aHr]± 


dr 

dvl 

dr 


d{aHr]j_ )^_^ 
dr 

d{aHr]±) rp 
dr 


+ O.TNVa + (HtT + k )Va = 0, 

-I- ^TNV^ + i^NN + k‘^)Va = 0- 


(3.18) 

(3.19) 
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We see that the 2-field perturbation system consists of a set of pairwise coupled harmonic 
oscillators. The coupling between curvature (uj) and isocurvature {v^) modes is controlled 
by r]± which is inversely proportional to the curvature radius of the background trajectory. It 
follows that the coupling between curvature and isocurvature will be strong whenever there 
is a sharp turn in the background trajectory. The symmetric mass matrix O of Eq. (3.18) 
and (3.19) has the following elements: 



+ 2e- 3r?|| -h ??||C|| - der/y 26^ - rj{) , 

(3.20) 

^NN = 

-a^H‘^{2 - e) + a^Vjvjv + a^H^eR , 

(3.21) 

^TN = 

a^H^r]±{3 + e-2r]\\ - C±) , 

(3.22) 


where R is the Ricci scalar of the scalar manifold and the third slow-roll parameters are 
defined as: 


'4>o 

U = (3.24) 

Assuming that the system is decoupled (r/x = 0) when observationally relevant modes 
are deep inside the horizon (^ 3> 1), canonical quantisation fixes the initial conditions for 
the scalar perturbations to be 


V 


N,T 

a. 




^—ikr 




(3.25) 


i.e. the initial conditions are simply given by the Bunch-Davies vacuum. Note that in Eq. 
(3.25) it is understood that = 1 etc. and derivatives define the initial condition for . 

It is important to note that there are two sets of equations of motion for the Muhkanov-Sasaki 
variables (in total 4 equations) and corresponding initial conditions, one for each value of a 
and that both should be taken into account when computing the inflationary observables. 
The power spectra are defined in terms of the scalar 2-point functions as 


^ , (3.26) 

^ a=l,2 

Ps{k,T) = v^{k,T)v^*{k,T) . (3.27) 


where P(^{k,T) and Ps{k,T) denote the dimensionless power spectra for the curvature and 
isocurvature modes respectively. Given that in multi-field models there can be superhorizon 
evolution of the perturbations, these are to be evaluated at the end of inflation. This is 
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to be contrasted to single field models, where the freezing of curvature perturbations on 
superhorizon scales means the power spectra can be evaluated at horizon exit. From Eq. 
(3.26) one can compute the spectral index for the curvature perturbations 


ris 


d ln{Pt^{k,Tend)) 
d ln(/i;) 


as well as the amplitude at the pivot scale fc* 


(3.28) 


■^s t'^end) t 


(3.29) 


which in the absence of an analytical solution to Eqs. (3.18), (3.19) must be computed nu¬ 
merically. Noting that the tensor modes’ evolution is unaffected by the number of dynamical 
fields driving the background expansion, the amplitude of the tensor power spectrum is given 
as in the single field case by 

At = Ptih,Tend) = , (3.30) 

TT 

which implies the following definition of the tensor to scalar ratio 


r = 


p< 


(fe*, Tend') 


2H^ 

TP(^{kif , Tend) 


(3.31) 


Besides probing the scalar and tensor power spectra, observations also put bounds on the 
total fraction of primordial isocurvature, defined as 


A, 


Ps 

Ps + Pq' 


(3.32) 


Erom the theoretical point of view the isocurvature fraction depends on the mass of the 
isocurvature modes {VLmn) and on the strength of their coupling to the adiabatic perturba¬ 
tions. The observational bounds on fdiso can vary by many orders of magnitude, depending 
on how primordial isocurvature is transferred to the post-inflationary Universe. From [58] 
we find that the less constraining bound is of the order 


Piso < 10-3 


(3.33) 


at the end of inflation. 

In addition to putting constraints on the fraction of primordial isocurvature, observa¬ 
tions also put constraints on the non-linear non-Gaussianity parameters Jnl- Producing 
large non-Gaussianities would, in principle, spoil the validility of the model. However, for 
2-field models the non-linear are of the order of the slow-roll parameters and hence they 
are surpressed in our model, see [59], [60]. We will not consider non-Gaussianities beyond 
this point. 
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3.1 Decoupling limit and the single field observables 


In order to understand how to relate the observables defined above with those of single field 
inflation let us once more take the decoupling limit ? 7 _l —)• 0 and reduce the equations of 
motion to 






+ 


+ 


k H—^ —2 + + 1 —3 + 


h H—2 (—2 — 06 + 3?7|| ) 
m2 f _ 2M2\ 

Ip j 


= 0 , 


= 0 , 


(3.34) 

(3.35) 


where the isocurvature mass is M^ = V)vAf + and we have used = aH{l — e) for the 

background evolution. Equations (3.34) and (3.35) admit solutions that are a superpositions 
of Hankel functions of the first and second kind, which upon imposing Bunch-Davies boundary 
conditions reduce to 


vj H^lH-kr) = v'^ , z/T = ^ + 2e - 77|| (3.36) 

V2 Hui^i-kr) ^ = (l + 2e) + ^e (3.37) 


up to overall unimportant phases, with uj = = 0. Note that vt = 3/2 and vn < 3/2 to 

zeroth order in the slow-roll expansion. 

On superhorizon scales, in the decoupling limit, one can show that the curvature per¬ 
turbations are frozen as in the pure single field case, since /a oc while the 

isocurvature perturbations decay as /a oc gyN-'ip ^ a~ to zeroth order in slow 

roll. 

In order to make contact with observations one can compute the dimensionless power 
spectrum of curvature perturbations in this limit, finding 


^ dvr^a^e \aHJ Svre 


(3.38) 


on superhorizon scales, which implies the following definitions for the spectral index and the 
amplitude 


Ug — 1 = 3 — 2 i't = —4 e + 2m = —2e — r/ and Ag = —^ . (3.39) 

" Svr^e 


Given that in the ri± = 0 limit the curvature perturbations are frozen once they leave the 
horizon, these observables can be evaluated at horizon exit, when k = aH. 

By performing a similar computation for the isocurvature modes one can show that 
the amplitude of the isocurvature power spectrum at horizon crossing is the same as for 
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the adiabatic modes. As noted above, due to the fact that vn < 3/2, the isocurvature 
perturbations decay on super horizon scales with a rate controlled by the ratio This 

implies that the isocurvature fraction at the end of inflation scales as 

Piso ~ TWTTo = a ^ (3.40) 

If-' p 

and is therefore suppressed if at some stage during observable inflation M > H. 

The tensor to scalar ratio is given in the decoupling limit by its single field expression 

r = ^ = 16e*, (3.41) 

-^S 

where * denotes evaluation at horizon crossing since both adiabatic and tensor perturbations 
are frozen outside the horizon. 

4 Results 

Having reviewed the formalism for the computation of inflationary observables in 2-field 
models we now reanalyse Higgs-otic inflation in three representative points in parameter 
space: A = 0.83, A = 0.7 and A = 0.2. First, we will keep the total amount of flux fixed at 
(5 = 1. For standard values of the string scale this leads to a supersymmetry breaking scale 
of Mss — 10^^ GeV. However, we will also show the results for a large range of G values at 
the end of this chapter. 

4.1 Higgs-otic regime: A = 0.83, (5 = 1 

In this section, we present the results for the inflationary observables corresponding to the 
canonical flux choice A = 0.83 and (5 = 1. For this point in parameter space, at the SUSY- 
breaking scale Mss — 10^^ GeV, there exist a heavy Higgs that can be integrated out and a 
light (approximately massless) Higgs, to be identified with the Standard Model Higgs field. 
As mentioned in Section 2, at the ultraviolet (string) scale the two helds have comparable 
masses (in fact for A = 0.83 one has mH/nih = 3.28), implying that inflationary dynamics 
driven by such a Higgs sector will necessarily be multi-field in nature and that the observables 
are better estimated via the methods reviewed in Section 3. 

We present in figure 3 several trajectories in the {a, 9) plane as well as the evolution of 
the all important r]± parameter for those same trajectories. Recalling that r]± is proportional 
to the inverse curvature radius of the background trajectory [55], we observe that the marked 
turns in the trajectories generate clear peaks in r]±. These peaks are sharper and higher for 
trajectories where the turning takes place close to the end of inflation (large 9o trajectories). 
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Figure 3. Left: The last 68 efoldings for some representative inflationary trajectories. Red dots 
mark 10 efoldings intervals on each trajectory. Right: Evolution of ri± for the different trajectories. 

The multi-field effects, arising through the ri± controlled coupling between curvature 
and isocurvature modes, are more severe the earlier the turn takes place. So even though the 
coupling is stronger for say trajectory C than for trajectory H, its effects are more pronounced 
in the latter since the isocurvature and curvature modes are coupled at earlier times, when 
there is more isocurvature power. 

The effects of the multi-field dynamics on the amplitude of the scalar curvature per¬ 
turbations are illustrated in figure 4 where we see that the power transfer from isocurvature 
to curvature is maximized the earlier the turn takes place. The impact of multi-field ef¬ 
fects on the scalar amplitude for trajectories A-E is minimal, with the single field estimate 
of Eq. (3.39) providing a good approximation to the full result. For trajectories F-H the 
superhorizon evolution of the curvature perturbations driven by isocurvature power transfer 
implies that Eq. (3.39) underestimates the amplitude by as much as 80%. These results are 
summarised in table 4.1, where we also present the estimates for the primordial isocurvature 
fraction, /3iso, at the end of inflation on the k^o and fego scales. We observe that f3iso varies by 
many orders of magnitude, being larger for late turning trajectories (large 9o) where power 
transfer between isocurvature and curvature perturbations is less efficient and the attenuation 
of isocurvature power is mostly driven by its decay on super horizon scales. 

In the same way as the multi-field effects can lead to an underestimate of the scalar 
amplitude, they will also impact other inflationary observables, in particular the tightly 
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Trajectory 

A 

B 

C 

D 

E 

F 

G 

H 


1.01 

1.01 

1.02 

1.03 

1.07 

1.22 

1.76 

1.716 

logic/?iso(A:6o) 

-3 

-3 

-5 

-8 

-13 

-16 

-18 

-20 


1.01 

1.01 

1.02 

1.04 

1.09 

1.30 

2.00 

1.33 

logic/3jso(A:5o) 

-3 

-3 

-5 

-8 

-13 

-16 

-18 

-19 


Table 1. Ratio between the amplitude of the curvature perturbations at the end of inflation and the 
single field estimate of eq. (3.39) for the trajectories of Fig. 3. 


constrained spectral index and the tensor to scalar ratio. The single and multi-field estimates 
for these quantities are plotted as functions of the initial condition 6q in figure 5. 

Starting with the tensor to scalar ratio, we observe that the effect of the multi-field 
dynamics is to flatten the peak and therefore to bring the results more in line with the 
PLANCK 2015 constraint of r < 0.114. This effect is partially due to the tensor modes 
beings unaffected by the multifield effects and in part due to the fact that the single field 
estimate for the amplitude of the scalar perturbations is a bad approximation for trajectories 
that turn early. Therefore by underestimating the amplitude of the scalar fluctuations, the 
single field formula overestimates the tensor to scalar ratio by 

r = 16 (A:*, Tend) ■ (4.1) 


In what concerns the spectral index, we also see that the multi-field estimate is con¬ 
siderably sharper than what one would expect by applying single field results. We observe 
that the peak at low Oq is absent and that the trough is shallower. We recall that the single 
field estimate for n^, Eq. (3.39), is obtained by taking the decoupling limit —)• 0. This 

condition is clearly violated by the early turning trajectories, for which ry_L peaks as the scalar 
modes leave the horizon. We therefore conclude that the peak in the low Oq region is spurious 
To understand the change in Ug for larger 9q trajectories it is useful to rewrite eq. (3.18) 
as 


'2„,r 






-b [VLtt - a -b k )v, 


7^ 

■^a 


elasticforce 


—2aHr]± 


K 

dr 


d{aHr]±) ^ 
dr 

■v* 

external force 


^TNVa ■ 


“^Comparing the results for the single field estimates of figure 5 and those of [18] we see that they differ in 

the low 9o range, where [18] has no peak. This difference can be traced back to how one generalises rj = ^ for 

V- • 

multifleld cases. If one takes rj to be the smallest eigenvalue of —then indeed there is no peak . However 
we use here a different, and more accurate, prescription which can be indeed derived from the decoupling 
limit, as we argued in section 3. 
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We see that the equation of motion for the curvature perturbation is equivalent to a fric¬ 
tionless harmonic oscillator with a ’’time” dependent proper frequency subject to an external 
force whose magnitude is set by the isocurvature perturbation. The effects of a turn in the 
background trajectory, which gives rise to the external force, will be more pronounced on 
/c-modes for which the ratio between the external force and the elastic force is larger. This 
ratio is well approximated by the simpler relation between 

m - fh . ( 4 . 2 ) 

Since the amplitude of different /o-modes around the pivot scale will be affected differently 
by a turn in the background trajectory, there will be super horizon evolution of the spectral 
index for the curvature perturbations. From the solutions in the decoupling limit {r]± = 0) 
one finds that on superhorizon scales and before the turn 


where denotes the curvature spectral index before the turn in the trajectory and is assumed 
to be < 1. As for the isocurvature perturbations, one may expand 

+ +t<^. (4.4) 

which leads to 

R{k) oc a > 0 . (4.5) 


To understand how this changes the spectral index consider a pair of k-modes around the 
pivot scale k^: k- < k^ < k+. It follows that since 


R{k-) 

R{k+) 


(X 



< 1 , 


the mode power is more enhanced than the k- mode power resulting in a spectral index 
closer to unity 


n 


end 

s 


> rVi 


(4.6) 


This behaviour can be clearly observed in figure 6, where the k- power is less enhanced than 
the A:+, resulting in a more even distribution of power and in an increase in the spectral 
index. Note that, before the turn (W < 50), the single field estimate is actually a good 
approximation to the full result and that it only fails due to the sharp turn in the background 
trajectory that causes conversion of isocurvature into curvature power. 
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Figure 4. Evolution of the curvature (blue) and isocurvature (red) two point functions for fceo for the 
trajectories of Fig. 3. Curvature and isocurvature power are normalised to the single field estimate 
P — I 

® S-TT^e I ’ 
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Figure 5. Spectral index (left) and tensor to scalar ratio (right) for (5 = 1 and A = 0.83 superimposed 
with the PLANCK 2015 constraints ris = 0.96688 zb 0.0061 and r < 0.114 (grey band). The slightly 
transparent curve corresponds to the single held estimate while the one surrounded by the data points 
corresponds to the multi-held results. Circles: = 60, Squares: Ne = 50. 




Figure 6. Time evolution of the curvature perturbations and of the spectral index on scales = fceo 

for the case A = 0.83,G = 1 and Oq = 0.9. 
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4.2 Modified Higg-otic regime: A = 0.7, G = 1 

The value of ^ = 0.83 analysed in the previons section follows, via RG running, from the as¬ 
sumption that the light scalar snrviving down to the EW scale is to be identified with the SM 
Higgs and that no extra degrees of freedom enter the particle spectrum beyond those of the 
MSSM. It is however conceivable that extra light particles are present in an extension of the 
MSSM, modifying the rnnning of the Higgses’ masses between the compactification/string 
scale and the SUSY-breaking scale M 55 . If this is the case then other valnes of A can be 
compatible with the Higgs-otic scenario. Keeping this possibility in mind we now analyse 
the case A = 0.7, which corresponds to a mass ratio mu/'fnh = 2.38 at the string scale. 




Figure 7. Left: The last 68 efoldings for some representative inflationary trajectories. Red dots 
mark 10 efoldings intervals on each trajectory. Right: Evolution of r]± for the different trajectories. 

In figure 7 we present sample background trajectories and the corresponding evolution 
of the r/_L parameter. Comparing with the results of the previous section we see that the 
trajectories are straighter and that the r]± peaks are less pronounced and located at later 
times. This implies that the differences between the exact results and the single field esti¬ 
mates for the observables should be less pronounced than for the A = 0.83 point. That is 
indeed the case, as can be seen from comparing figures 5 and 8 . Though less pronounced, the 
disparity between single field estimates and the full results is still important as the effect of 
the multi-field dynamics is to bring the observables more in line with the current constraints 
on Hs and r: the variation in the spectral index is damped and the tensor to scalar ratio is 
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significantly reduced. In fact Ug is comfortably inside the 2(t PLANCK 2015 band and the 
r is reduced to the point that it complies with the observational upper bound V 6q. This 
highlights the importance of properly estimating the observables at a time of ever increasing 
measurement precision. 
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Figure 8 . Spectral index (left) and tensor to scalar ratio (right) for G =1 and A = 0.7 superimposed 
with the PLANCK 2015 constraints Ug = 0.96688 ± 0.0061 and r < 0.114 (grey band). The slightly 
transparent curve corresponds to the single field estimate while the one surrounded by the data points 
corresponds to the multi-field results. Circles: Nf. = 60, Squares: Nf. = 50. 

In table 4.2 we present the isocurvature fraction at the end of inflation as well as the 
comparison between the single field estimate for the adiabatic amplitude and the numerical 
result for the various trajectories. 

4.3 Almost single field regime; A = 0.2, G = 1 

The Higgs-otic idea also features regimes in which the connection between inflation and 
MSSM Higgs physics is absent or hard to realise. In these cases, where the fluxes are such 
that A is substantially different from the canonical value of 0.83, inflation can still be driven 
by the D-brane position modulus on the compact space’s if the inflaton is associated to 
other degrees of freedom. Indeed, for small A the values of the Higgs masses m/i, mn are too 
close to each other for the running from the string scale down to the SUSY-breaking scale 
to be sufficiently strong to yield an (approximately) massless SM doublet at M 55 , but one 
could still identify the complex inflaton with other degrees of freedom in some extension of the 
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Trajectory 

A 

B 

C 

D 

E 

E 

G 

H 


1.01 

1.01 

1.02 

1.04 

1.09 

1.20 

1.40 

1.22 

logio/3iso(/C6o) 

-4 

-4 

-5 

-9 

-10 

-9 

-12 

-13 


1.01 

1.01 

1.02 

1.05 

1.11 

1.26 

1.43 

1.16 

logio Pisoiho) 

-3 

-3 

-5 

-9 

-9 

-11 

-12 

-13 


Table 2. Order of magnitude of the isocurvature fraction at the end of inflation for the trajectories 
of Fig. 7. 

MSSM. Note, in particular, that if SUSY particles are found at LHC, the canonical Higgs- 
otic scenario with A = 0.83 would be ruled out, since it assumes a large SUSY-breaking 
scale M 55 ~ 10^^ GeV. In this case, the inflaton could however be identified with extra 
scalars which could have SUSY preserving masses at the inflaton scale ~ 10^^ GeV. As we 
said, examples of those degrees of freedom could be right-handed doublets of a left-right 
symmetric extension of the MSSM or even colour triplet scalar partners of the Higgs fields. 
In such a case, very small values for the inflaton parameter A would be consistent. 

While the search for specific MSSM extensions with this structure is quite interesting, 
the inflationary dynamics may be studied in a model independent manner assuming that 
such new non-Higgs degrees of freedom correspond again to the position moduli of D-branes. 
We study here for comparison the case with A = 0.2 (corresponding to mn/Trih = 1.22). 

In figure 9 we present the background evolution for sample trajectories in such a regime. 
In this case we see that the trajectories are essentially straight along the a direction with 
only slight turning in the last 10 efoldings of expansion. This is in accordance with the fact 
that 9 becomes massless in the limit of vanishing A. The straight trajectories imply that t]±, 
being inversely proportional to the curvature radius, vanishes everywhere except at the very 
end of inflation (where the mild turning takes place) as can be seen in figure 9. As mentioned 
above rj± controls the coupling between the curvature and isocurvature modes and so this 
limit corresponds to a decoupling regime of Eqs. (3.18) and (3.19). We therefore expect the 
single field estimates of Eqs. (3.39) and (3.41) for the inflationary observables to provide a 
good approximation to the full result. In fact, if one employs the multi-field formalism in 
the computation of the observables and compares it with the single field estimates, one finds 
that there is agreement at the level of a few percent. This is displayed in figure 10 where the 
naive single field bands track the exact results to the point of being almost indistinguishable. 

Since isocurvature is practically decoupled from curvature throughout the observable 
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Figure 9. Left: The last 68 efoldings for some representative inflationary trajectories. Red dots 
mark 10 efoldings intervals on each trajectory. Right: Evolution of ri± for the different trajectories. 

inflationary range, the isocurvature fraction at the end of inflation as estimated by /3jso is 
larger that in the previous cases, as can be seen in table 4.3. This is due to the fact that, 
unable to transfer power to the adiabatic mode, all the isocurvature modes can do in their 
superhorizon evolution is to decay slowly. Note that even though this almost single field 
regime gives rise to the largest isocurvature fraction, it is still below the highest upper bound 
on Piso derived from the latest PLANCK 2015 data [58]. 


Trajectory 

A 

B 

C 

D 

E 

F 

G 

H 


1.01 

1.01 

1.01 

1.02 

1.02 

1.02 

1.02 

1.01 

^Ogio Pisoiho) 

-4 

-4 

-5 

-4 

-5 

-6 

-6 

-6 


1.01 

1.01 

1.02 

1.02 

1.03 

1.03 

1.02 

1.02 

logio/3iso(/c5o) 

-4 

-4 

-4 

-4 

-5 

-5 

-6 

-6 


Table 3. Ratio between the amplitude of the curvature perturbations at the end of inflation and the 
single field estimate of eq. (3.39) and amplitude of isocurvature perturbations for the trajectories of 

fig- 9- 


- 25 - 



























0.99 

0.98 

0.97 

=” 0.96 

0.95 

0.94 

0.93 


0.0 0.5 1.0 1.5 2.0 2.5 3.0 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

00 


0.99 

0.98 

0.97 

0.96 

0.95 

0.94 

0.93 


0.0 0.5 1.0 1.5 2.0 2.5 3.0 

0.20 

0.18 

0.16 

0.14 

0.12 

0.10 

0.08 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 


0.18 


0.16 


0.14 


0.12 


0.10 - 


0.08 


0(1 


Figure 10. Spectral index (left) and tensor to scalar ratio (right) for G = 1 and A = 0.2 superimposed 
with the PLANCK 2015 constraints ng = 0.96688 ± 0.0061 and r < 0.114 (grey band). The slightly 
transparent curve corresponds to the single field estimate while the one surrounded by the data points 
corresponds to the multi-field results. Circles: Ng = 60, Squares: Ng = 50. 


4.4 Varying G 

Finally let us discuss the effects of varying G over the results. As we commented in section 
2, G parametrises the total amount of flux quanta in Planck units, which determines in turn 
the ratio between the supersymmetry breaking scale and the string scale. In particular, 


G = ~ 5 

h 


Mss 

Mr, '' M2 


(4.7) 


with n being the flux quanta. We will take the string scale of order Mg r\j 10^® - 10^^ 
GeV, as suggested by gauge unification. An isotropic compactification with n ~ 0{1) then 
implies Mss ~ 10^^ GeV. But it should be noticed that the parameter n in general receives 
contributions from a large number of 3-cycles so that large cancellations can take place 
that lead to n <C 1, lowering the scale of SUSY breaking. However, n <C 1 is problematic 
if one wants to satisfy the experimental bounds on density scalar perturbations. A lower 
supersymmetry breaking scale may lead to a too low amplitude of the scalar power spectra 
in this inflationary model. The best fit indeed corresponds to Mss ~ 10^^, which is the 
typical scale for flux-induced supersymmetry breaking obtained by assuming the flux quanta 
to be of order one. On the other hand, n S> 1 (and thereby a higher supersymmetry breaking 
scale) makes the potential energy trans-Planckian at the beginning of inflation, which is 
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inconsistent with the effective field theory approach. Throughout this paper we have shown 
the results for the representative value (5 = 1 in Planck units. Here we also show the results 
for other possible values of G for completeness. In figure 11 we plot our results in the r — Ug 
plane superimposed over the experimental Planck exclusion limits. The data correspond to 
A = 0.83 and arbitrary G. 


s 



0.25 


Primordial tilt 


Figure 11. Tensor to scalar ratio and spectral index for Higgs-otic inflation with A = 0.83 and 
arbitrary G. The data is superimposed over the recent Planck exclusion limits [58]. The color pattern 
(from red to blue) corresponds to (higher or lower) density of initial condition points . There is a 
lower cutoff on the density required to be plotted (or equivalently in the level of fine-tuning allowed) 
missing around 10 % of the points. 


Notice that G not only enters in the absolute value of the potential, but also in the field 
redefinition to get canonical kinetic terms. The bigger G is, the stronger is the flattening 
of the potential and the results are closer to those of linear inflation. This is the reason 
why our results actually interpolate between quadratic (small G, negligible flattening) and 
linear (big G, strong flattening). However, as we already said, those closer to linear inflation 
are in better agreement with both density scalar perturbations and tensor-to-scalar ratio 
constraints from Planck. In figure 11 the color pattern from red to blue refers to the density 
of points, being the red regions the most populated. This could have been anticipated from 
figure 5, where most of the initial conditions gave rise to values of r and closer to the single 
field prediction. The spreading of the results to smaller values of the spectral index is due to 
the freedom on the choice of the initial conditions, but as we said the blue region corresponds 
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to very fine-tuned values of 9q and the majority of the points is localised at Ug ^ 0.965 and 
r ~ 0.08 — 0.12. It is remarkable that the most populated regions are indeed those in better 
agreement with experiments. 

5 Conclusions 

The current lack of experimental evidence for low-energy supersymmetry coupled with the 
possibility of observing primordial gravitational waves in the near future drives us to analyse 
a model in which these two seemingly unconnected phenomena have their roots in the same 
sector of the theory: Higgs-otic inflation. 

In Higgs-otic inflation the scale of SUSY breaking, Mss ~ 10^^ GeV, is identified with 
the inflation’s mass and the Universe’s expansion is driven by the MSSM Higgs scalars. Such 
a high inflationary scale gives rise to a sizeable tensor fraction and implies super-Planckian 
field excursions for the inflatons. The model is embedded into a local string construction, in 
which the Higgses/inflatons are D-brane position moduli. Such a stringy embedding allows for 
control over higher order corrections to the inflationary potential that are usually problematic 
for large field inflation models. 

In this work we have analysed in detail the observational signatures of this model, taking 
into account the effects of isocurvature perturbations that are unavoidable whenever there are 
multiple dynamical fields during inflation. We found that the effects of 2-field dynamics on 
the observables are more pronounced whenever there is a turn in the background trajectory. 
In Higgs-otic inflation the presence of such a turn during the last 60 efoldings of expansion 
is controlled by the choice of background fluxes. For flux choices leading to almost straight 
trajectories, the single field estimates for the inflationary observables derived in [18] are a good 
approximation to the full result as we have seen in Sec. 4.3. On the other hand for those flux 
choices consistent with the identification of the inflaton with the MSSM Higgses, we find that 
the trajectories will necessarily present some degree of turning. At the turn isocurvature and 
curvature perturbations are coupled, with isocurvature sourcing the curvature power. This 
implies superhorizon evolution of the curvature perturbations and translates into an increase 
of the adiabatic perturbations’ amplitude and the corresponding decrease in the tensor to 
scalar ratio. Since different fc-modes are affected differently by a turn in the trajectory, there 
will be superhorizon evolution of the spectral index, which we have shown to lead to a sharper 
prediction for Ug when compared to the single field approximation employed in [18]. In what 
concerns the isocurvature power at the end of inflation, we have shown that it is always 
suppressed enough to be compatible with the less stringent upper bounds from PLANCK 
2015. We have shown in 4.1 and 4.2 that the observational footprint of Higgs-otic inflation 
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is narrower than what one would expect by simply performing an adiabatic projection: the 
range over which Ug varies is smaller and centred around the Planck best fit value and the 
tensor to scalar ratio is reduced significantly yielding a region r = 0.08 — 0.12 depending on 
the model initial conditions. 

In summary, our analysis demonstrates that the observable signatures of Higgs-otic 
inflation are in line with the results of the joint Planck/BICEP analysis [61]. It would be 
interesting to address additional questions like moduli stabilisation and reheating in this 
model. We leave those aspects to future research. 
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